Universal relation between the dispersion curve and the ground-state correlation 
length in ID antiferromagnetic quantum spin systems 
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We discuss an universal relation s(iKj = with Ken = l/£ in ID quantum spin systems with 
an excitation gap, where e(k) is the dispersion curve of the low-energy excitation and £ is the 
correlation length of the ground-state. We first discuss this relation for integrable models such as 
the Ising model in a transverse filed and the XYZ model. We secondly make a derivation of the 
relation for general cases, in connection with the equilibrium crystal shape in the corresponding 2D 
classical system. We finally verify the relation for the S — 1 bilinear-biquadratic spin chain and 
5 = 1/2 zigzag spin ladder numerically. 



I. INTRODUCTION 

Recently, considerable attention has been paid for the 
low-energy excitations of one-dimensional(lD) antiferro- 
magnetic(AF) quantum spin chains with an excitation 
gap. In discussing the low-energy excitations of the spin 
gap systems, one of the most fundamental quantities is 
the dispersion curve e(k) of the elementary excitation. 
We can then expect that e(k) reflects the nature of the 
ground state correlation function G(n) = (SfiS*), since 
the elementary excitation can be interpreted as a "lo- 
cal defect particle" in the ground-state. Actually for the 
5=1 AF Heisenberg chain, a relation v — £A is known 
numerically, where v is the spin wave velocity, £, is the 
correlation length and A is the excitation gapld How- 
ever, a full justification of such a relation has not been 
made yet and to find its generalized form is an interesting 
problem. 

In this paper, on the basis of the correspondence 
between the ID-quantum and 2D-classical systems, we 
present an universal relation between the low-energy dis- 
persion curve and the ground-state correlation length: 



e{iK) = with Rck = l/£, 



(1) 



which holds for a wide class of ID quantum spin chains. 
A remarkable point in eq. (pi), which has not been 
stressed so far, is that the full dispersion curve should 
be taken into account properly. In other words, the con- 
ventional relation v = £A does not hold for general cases; 
v — £A is valid only for the relativistic-free-fermion case, 
where the dispersion curve is e(k) = y A 2 + v 2 k 2 with 
k = it — k. Although the form of eq.(0) looks quite nat- 
ural on the assumption of the free quasi-particle theory, 
the validity of such a free particle picture is a highly non- 
trivial problem in the strongly correlated systems. 

In the following, we address the relation (|l|) both in 
analytical and numerical ways. We firstly show that the 



relation is actually satisfied for integrable models, where 
we can handle the ID-quantum and 2D-classical corre- 
spondence exactly. This analysis of the integrable mod- 
els provides an essential insight for general cases. We 
secondly make a derivation of the relation for the gen- 
eral cases, dealing with the Suzuki- Trotter decomposed 
2D classical systemsu. By using density-matrix renor- 
malization group(DMRG)B and the exact diagonalization 
tcquniches, we finally demonstrate that the relation is ac- 
cepted for the 5=1 bilinear-biquadratic(BLBQ) chain 
and the S = 1/2 zigzag spin ladder, which exhibit char- 
acteristic dispersion curves associated-with the cusp sin- 
gularity in the magnetization process.ETEl 



II. EXACTLY SOLVABLE CASES 

Although a lot of important relations are established 
in the integrable systemsEI, eq. (|l|) has not been dis- 
cussed yet. Thus we start with analyzing some inte- 
grable lattice models. As the simplest but important 
example, let us examine the ID Ising model in a trans- 
verse field(ITF): H = JJ2 s i s i+i ~ ^J2 s fy where s t is 
the spin- 1/2 operator at i th site and F is the trans- 
verse field. This model is solved with the Jordan- Wigner 
transformation, and then the dispersion curve is o. 



tained as £i s in g (fc) = \\J{J cosk + 2T) 2 + (Jsinfc) 2 .t3 
The correlation function can be also calculated straight 
forwardly to be £ -1 = |ln(2J/r)|. We can then confirm 
£isin g (J/C) = clearly. 

More generally, we consider the XYZ model: H = 
J2i JxO-fa-f +1 +J y a^a^ +1 +J z afal +1 , where <7j is the Pauli 
matrix at i-th site and the coupling constants are repre- 
sented by the Baxter's elliptic- function parametalization 
with a modulujB I; J x = cn2£/sn2£, J y — dn2C/sn2C and 
J z = l/sn2£.0 We consider only the AF region of the 
couplings. The dispersion curve of the spinon excitatioe. 
is calculated exactly by Johnson, Krinsky and McCoy:liil 



£xYz(fc) = 2 J; 



K lS n{2{) 
K' 



1 — k\ cos 2 k, 



(2) 



where K' is the complete elliptic integral of the modulus 
I, and fci is the modulus defined by K[/Ki = irQ/K' . The 
correlation function of the XYZ chain is also calculated: 



G{n) 



d4>id<t>2p(4>i,4>2) 



K 



«2Sn( )sn( 



K 
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(3) 



where p(4>i,4'2) is a non-singular function and fe is an- 
other modulus defined by the elliptic integrals K^fKi = 
2-kC,/K' Jill In order to compute the asymptotic form 
of the correlation function, we apply the saddle point 
method to eq. (g) and obtain G(n) ~ fc^. Since the a z 
operator creates two spinon excitations when it acts on 
the ground state wavefunction, the spinon-spinon corre- 
lation length is given by £ _1 = — ^lnfo. Then, we can 
see £xyz(^- lnfe) = easily, by using the Landen trans- 
formation. From these, we have verified eq. (El) for the 
integrable models. 



III. GENERAL CASES 

The realization of eq. (0) for the integrable models can 
be attributed to the nature of the underlying 2D classical 
models, such as 8-vertex model. An useful point in con- 
sidering the 2D classical systems is that the asymptotic 
behavior of the system in the time and space directions 
can be treated systematically; Particularly, it should be 
remarked that the anisotropic correlation length of the 
system can be associated with the equilibrium-crystal 
shape(ECS) of the interface via the Wulff construction:E3 



7(0) = 



k B T 

£(0)' 



(4) 



where 7(6*) is the crystal shape at the angle 6 = 
arctan(y/x) and £((9) is the angle dependent correlation 
length. For the integrable models, the universal behav- 
ior of the ECS is shown exactly, and then the eq. (|l|) for 
the corresponding ID quantum spin systems can be ex- 
plained as one of the resulting behaviors of the universal 
ECS. 

In order to see the connection of the ECS with the 
quasi-particle excitation, we represent the correlation 
length £(0) as the asymptotics of the lattice Green's func- 
tion: 



— — = - lim ~1nG(x,y) 

£(p) r^oo r 



(5) 



where G(x,y) is the lattice greens function and r = 
y ' x 2 + y 2 . Throught the Fourier space pepresentation, 
we can write this lattice Green's function as 



G(x,y) 









JJ(k x , Ky) 



(6) 



in which the asymptotic behavior in the reaLspace is de- 
termined by the condition D(k x ,k y ) — O.Ejo The corre- 
lation length and the dispersion curve in the ID quantum 
system can be also extracted from D(k x ,k y ) (see below 
in detail). In_fact D(k x ,k y ) of the XYZ model is calcu- 
lated in Ref.Ej, where we can reduce the universal curve 
of the ECS toeq. (Q). 

In the framework of the Wullf construction, what we 
want to emphasize here is that there is the one-to-one 
correspondence between the ECS and the asymptotics of 
the lattice Green's function. Since the ECS is the well- 
defined thermodynamic objecttS and the Wulff construc- 
tion does not require the integrability of the system, we 
can assume the existence of the lattice Green's function 
with D(k x ,k y ) for a wide class of non-integrable mod- 
els. This implies that the asymptotic form of the lattice 
Green's function (|6|) can be defined not by the perturba- 
tive argument of the local excitation, but by the thermo- 
dynamic quantity of the ECS, where the interaction effect 
of the system is fully taken into account. 

To analyze the general ID quantum systems, we 
consider the dynamical correlation function G{n, r) = 
(S z (n, t)S z (0,0))i3, where S z (n,r) is the spin operator 
of the Heisenberg representation at position n and imag- 
inary time r, and (• • -)/3 means the thermal average at 
an inverse temperature j3. Our goal is to evaluate the 
asymptotic behavior of G(n, r) in the zero-temperature 
limit. For this purpose we describe G(n, r) by the lattice 
Green's function dq) on the Suzuki- Trotter transformed 
2D classical system of the size N x 2M: 



G(n, 2m) = / dk x 






a ik x n-\-ik y 2m 



D\K X , Ky) 



(7) 



where m = Mt//3, N is the lattice size in the spatial 
direction, and M is the Trotter number. In the Trotter 
limit M — * 00, we recover the original Green's function 
G(n,r) ~ G(n,2m). The function D(k x ,k y ) should be 
determined when a model is given. For example, that 
of the ITF is given by D(k x , k y ) — cosh 2^ cosh 2K y — 
sinh 2K X cos k x — sinh 2K^, cos k y with K x = (3 J /(AM) 
and Ky — i ln[coth(^)].E 4 J However, the explicit form of 
D(k x , ky) is not required in the following discussion. 

In eq. M), we first consider the Trotter limit M —* 
00, where the dominant contribution to the k y integral 



comes from k y ~ 0, since e 



ik v 2m 



is highly oscillating. 



Substituting k y — ^fc s , we can expand D(k x ,k y ) as 



a 



D{k X ,ky) - D(k X ,0) + Dyy{k X ,Q){^kyY + ■ ■ ■ , (8) 

where the first-order term of the k y derivative vanishes 
due to the Helmiticity of the original ID spin chain. 
Though D{k x ,k y ) contains the Trotter number depen- 
dence, we find that D(k x ,0) + Dyy{k x ,0)(j^k y ) 2 has a 
proper limit as M — > oo, without the overall normaliza- 
tion factor, and then the higher-order terms with respect 
to k y can be neglected. Thus we obtain 



g(n,2m)= / dk x F(k x ) 



(In,, 



ik x n-\-ik v T 

' d{k x ) + k~l 



(9) 



where F{k x ) is a non-singular function and we have de- 
fined 



d(k x ) 



2MD(k x ,0) 
m-^oo j3D yy (k x ,Q) ' 



lim 



(10) 



For the case of the ITF, we have d(k x ) = \{J 2 + AT 2 + 
4JT cos k x ), and F(k x ) = 2MT//3 explicitly. 

We next consider the low-energy limit (3 — ► oo. We 
then evaluate the "long-distance behavior" of eq. (Q) 
with a fixed direction r/n = tan# in the limit n, r — > oo. 
Applying the saddle point method, we obtain G(n, t) ~ 
where Kj and K a satisfies the saddle point 



&X ** fcyT 



conditions: 



d{iK x ) 



2k, 






tan 6*. 



(11) 



Particularly for the case of the equal-time correlation 
function, we set Kj, = 0(0 = 0). Then G{n,r) reduces 
into G(n) ~ e _Ka; ™ with d(in x ) = 0, in which the real 
part of k x represents the inverse correlation length of the 
ground-state. 

Let us here discuss the meaning of the function d{k x ). 
To see it, we further perform the k y integral in eq. (H) and 
then obtain another expression of the Green's function: 



G(n,r) 



const 









-yfdikx 



\Jd{k x ) 



(12) 



On the other hand, the asymptotic behavior of the 
Green's function can be reduced into G(n, r) ~ 
j dkG(k)e tkn ~ £ ( k " >T in the low-energy limit, by taking ac- 
count of the one-particle part of the excitation, where 
G(k) is the Fourier transform of G(n). Comparing the 
above expression with eq. (]12|), we can see \Jd{k) — e(k), 
which is checked for the ITF easily. Hence, we finally ar- 
rive at the formula M). 



IV. NUMERICAL RESULTS 

In the above derivation, we have formally written down 
the function D(k x , k y ). In this section, we verify the re- 
lation (|l|) for non-integrable models in numerical ways, 
since the analytic calculation of D{k x ,k y ) is very diffi- 
cult. As a typical example of the non-integrable models, 
we investigate the S = 1 BLBQ chain, which plays an 
important role to understand the Haldane systems. The 
Hamiltonian of the BLBQ chain is given by 



H = 



n 



Si ■ Si+i + (3(Si ■ Si+if 



(13) 



where Si is the 5 = 1 spin operator at i-th site, and (3 
is the coupling constant of the biquadratic term(not in- 
verse temperature). The BLBQ chain with (3 — 0, 1/3 
and 1 corresponds to the S — 1 Hcisenberg spin chain, 
the valence-batLd-solid(VBS) chairO and the SU(3) in- 
tegrable chaiixJ respectively. 

In connection with the relation (pi), a remarkable fea- 
ture of the BLBQ chain is that its dispersion curve 
and the ground-state correlation function exhibit various 
characteristic properties; Some extensive numerical stud- 
ies of the ground-state correlation function claim that 
the VBS point is a disorder point, which is the onset 
of the incommensurate oscillation-pf G{n) accompany- 
ing with a cusp singularity in £.EJ Appearance of Lif- 
sitz point (3 = (3l — 0.43 is also reported, where the 
peak wavenumber of-the static structure factor starts to 
shift as /? increases Jl3 Moreover, at (3 = (3 C ~ 0.40, the 
shape of the dispersion curve changes into the double well 
structure Jj which is associated, with the cusp singularity 
in the magnetization process B Although a phenomeno- 
logical explanation of this scenario is i«<;ently made by 
assuming a continuous-free-field theoryj23 a full justifica- 
tion based on the lattice Hamiltonian ([13]) has not been 
achieved yet. We now demonstrate the our relation (fil) 
for the BLBQ chain, which provides the background for 
the phenomenological theory. 

The ground-state correlation function of. the BLBQ 
chain can be investigated with the DMRG,E£l while there 
is few systematic study of the whole shape of the disper- 
sion curve for the sufficiently long chain. For the pur- 
pose of computing the dispersion curve, we employ the 
DMRG combined with the continued fraction method, 
where the excitation spectrum is extracted as a peak of 
the imaginary part of the dynamica l c or relation function 
expressed as the continued fraction.EjEa The DMRG cal- 
culation is performed for the BLBQ chain of the length 
L = 240 with the open boundary condition, targeting 
up to the five Lanczos vectors generated form the ini- 
tial state S^O), where £| = £S*e ifcn . In addition to 
this, we use the filtering technique to cut-pff the scatter- 
ing at the boundaries, according to Refo. In FIG.l we 
show the dispersion curves for (3 — 0, 1/3 and 0.6, which 
are calculated with the maximum number of the retained 
bases m = 81. Clearly we can see the shape change of 
the dispersion curve induced by the incommensurate ef- 
fect originating from the biquadratic term. The dotted 
parts in the curves for (3 = 1/3 and 0.6 mean that the dis- 
persion curves are screened by the continuum spectrum. 




FIG. 1. Dispersion curves of the BLBQ chain for /3 = 0, 
1/3, and 0.6 

In order to check eq. (||) , we estimate the Fourier series 
expansion of the dispersion curve with a fitting. Using 
the fitted curve, we determine the zero of the dispersion 
curve in the complex plane and compare the result with 
the correlation length obtained from a direct computa- 
tion of G(n). Here let us recall the dispersion curve has 
the square-root form e(k) = ^J 'd(k). Thus we can expect 
that the function d(k) is fitted well rather than e(k): 



N 



d(k) = e(k) 2 = 2_) A n cos nk, 



(14) 



71=0 



where {^4 n } is fitting parameters and N denotes the 
number of the parameters(cut off frequency). For the 
Heisenberg point (3 — 0, the fittings with iV_j= 5, 6 
and 7 yield £ = 0.590(5), 6.02(8) and 0.60(1) J_3 which 
are converged with respect to N and in good agreement 
with the known value £ = 6.03. We also test the fit- 
ting: s(k) = B + Bi cos(fc) + B 2 cos 2k H , and obtain 

£ = 2.5(1), 2.8(2) and 3.2(2) for N = 6, 7 and 8 re- 
spectively, which does not converge with respect to N. 
This slow convergence for e{k) may be essential, since 
e{Ik) — corresponds to the blanch-cut singularity of 
the square-root function. Hence, we have verified the 
e(iK) = with k = l/£ — in, including the square-root 
form of the dispersion curve. Here, it should be noted 
that we obtain £ = 6.00 for the case of assuming the rela- 
tivistic free fermion: e(k) — \J A 2 + v 2 k 2 with A = 0.410 
and v = 2.46. 

At the VBS point, the fittings with N = 3, 4 and 5 
yield £ = 0.92(1), 0.90(1) and 0.91(1) respectively, which 
are consistent with the exact value £ = 1/ In 3 = 0.910 
On the other hand, the relation v = £A with A = 0.70 
and £ = 0.91 gives the velocity v — 0.64, which does 
not agree with the correct value v = 0.85 estimated 
by the curvature at the bottom of the dispersion curve. 
Thus the conventional relation is not maintained in gen- 
eral cases. For the VBS model, we can further dis- 
cuss e(n + i/£) = exactly, based on the domain-wall 
picture of the excitation, where the dispersion curve is 
calculated as the Fa-acier transform of the domain-wall 
correlation function.cilE3 Since the asymptotic behavior 



of the domain-wall correlation function is determined 
by the overlap integral between the four degenerating 
VBS states, we can see the dispersion curve contains 
the factor 5 + 3 cos A:, which yields the zero of e(k) at 
k = n + z In 3. In addition, we should remark that the 
solution of d(k) — is given by the twofold root of 
k = n + iln3, because d(k) = e(k) 2 . This implies that 
the VBS model locates at the disorder point, since the 
imaginary part of the solution of d(k) = starts to move 
from 7T with making a complex-conjugate pair, as (3 is 
increased. 

For (3 = 0.6, the numerical data in the range ir/2 < 
k < n is well converged, where the imaginary part of the 
dynamical correlation function has sharp peaks. On the 
other hand, the curve in 0.67T < k < n/2 is disturbed 
by the continuum spectrum. However, we find that the 
main contribution to the zero of the e(k) comes from the 
range tt/2 < k < n. Actually, we have obtained the suffi- 
ciently stable results: £ = 3.79(5) and q — 0.678(l)7r and 
£ = 3.80(2) and q = 0.680(2)tt, by the fitting of eq. @) 
with N = 7 and 8 respectively, where q denotes the in- 
commensurate period of the correlation function. These 
values agree with £ = 3.87 and q = 0.6787T, which is com- 
puted directly from the real-space correlation function. 
Thus the formula (0) is verified for the BLBQ chain, in- 
cluding the period of the incommensurate oscillation. 

The similar behaviors of the dispersion curve and the 



in a class of the 
3 In the remaining 



correlation function can be expet 

frustrated spin ladder systemsElH! 

part of this section, we further test the relation (Q) for the 

5 = 1/2 zigzag spin ladder in the dimer gap phase, which 

also shows the characteristic behavioES of the dispersion 

curve and the correlation function.Q'E3 The Hamiltonian 

of the zigzag ladder is given by 



H 






Si ■ Si+i + aSi ■ Si+2 



(15) 



where S is the S = 1/2 spin operator. The zigzag ladder 
is in the dimer gap phase for a > 0.2411, and a = 0.5 
is called Majumdhar-Gosh point, where the ground state 
is exactly represented as the dimer singlet state. For 
a > 0.54, the spinon dispersion curve has the double 
well structure, in accordance with the incomennsurate 
oscillation of the correlation function. The magnetizati 
curve has also a cusp singularity in this region of a£rE3 
An important feature of the zigzag ladder is that the 
elementary excitation is described by the topological ex- 
citation of S — 1/2, which is called spinon. Recently, the 
single spinon excitation of the zigzag ladder has been ex- 
tracted by using the exact diagqiialization method with 
the Mobius boundary condition^ (See Fig. 2) although 
it has not been obtained for the system with the periodic 
boundary condition where the excitation is basically de- 
scribed by the combination of the two spinons. From 
the spinon dispersion curve calculated with the Mobius 
boundary condition, we estimate the square of the disper- 
sion curve d(k) with the fitting of (|l4|) and solve d(in) = 



numerically. We then compare the results with the cor- 
relation length computed with the DMRG method. 
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FIG. 2. Single spinon dispersion curves of the zigzag spin 
ladder for a = 0.48(open circle) and 0.6(solid circle). 

We here consider the zigzag ladder of a — 0.48 and 
a = 0.6 typically, whose properties are similar to the 
BLBQ chain of /3 = and (3 — 0.6 respectively. For 
a = 0.48, the period of the correlation function is com- 
mensurate and thus the bottom of the dispersion curve 
is located at k = tt/2. In fact, the result of d{in) = 
yields £ = 5.71(6), which is consistent with the correla- 
tion length obtained with the DMRG: £ = 5.68(7). Fol- 
ey = 0.6, the single spinon dispersion curve has a double 
well structure around k = tt/2, and then the solution of 
d(in) = yields £ = 4.35(20) and q = 0.368(5)tt. On the 
other hand, the asymptotic behavior of the correlation 
function computed with the DMRG is £ = 4.26(20) and 
q = 0.367(2), which is also in good agreement with the 
results of the relation |l]). We have therefore verified the 
relation (pi) for the S = 1/2 zigzag spin ladder, includ- 
ing the period of the incommensurate oscillation of the 
correlation function. 



V. SUMMARY AND DISCUSSION 

We have presented the relation ([!]) , which connects the 
zero of the dispersion curve with the correlation length 
of the ground-state. For the integrable systems, we have 
verified the relation by the exact expressions of the dis- 
persion curve and the correlation length. We have next 
derived the relation (|l|) for general cases, based on the 
Suzuki- Trotter transformed 2D classical models, where 
the connection to the universal crystal shape of the equi- 
librium interface plays an important role. To check eq. 
(0) in the non-integrable systems, we have computed the 
dispersion curve of the BLBQ chain, using the DMRG 
method combined with the continued fraction method. 
We have also investigated eq. (|l|) for the 5=1/2 zigzag 
spin ladder in the dimer gapped phase. Then we have 
shown that the relation (Jy) explains the asymptotic be- 



havior of the correlation function including the incom- 
mensurate oscillation. 

Since the argument in this paper is founded on the 
equilibrium crystal shape in the corresponding classical 
system, the validity of the relation (p]) is independent of 
such special properties of the system as the integrabil- 
ity, the anisotropy of the system(full-S't/(2) symmetry), 
the spin magnitude of the excitation particle, and so on. 
Indeed, the qasi-particle excitations considered in this 
paper includes various types of the excitations. We can, 
therefore, conclude that the relation (pi) is hold for a wide 
class of quasi-lD quantum systems. Recently, some inter- 
esting problems between the correlation function and the 
dispersion curve are also expected in a class-qf frustrated 
spin systems or correlated electron systems EQ When ana- 
lyzing such systems, we believe that the relation (nl) plays 
an important role. 
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